Phenotypic Characterization of Recombinant Marek's Disease Virus in Live Birds Validates Polymorphisms Associated with Virulence.

Marek's disease (MD) is a highly infectious lymphoproliferative disease in chickens with a significant economic impact. Mardivirus gallidalpha 2, also known as Marek's disease virus (MDV), is the causative pathogen and has been categorized based on its virulence rank into four pathotypes: mild (m), virulent (v), very virulent (vv), and very virulent plus (vv+). A prior comparative genomics study suggested that several single-nucleotide polymorphisms (SNPs) and genes in the MDV genome are associated with virulence, including nonsynonymous (ns) SNPs in eight open reading frames (ORF): UL22, UL36, UL37, UL41, UL43, R-LORF8, R-LORF7, and ICP4. To validate the contribution of these nsSNPs to virulence, the vv+MDV strain 686 genome was modified by replacing nucleotides with those observed in the vMDV strains. Pathogenicity studies indicated that these substitutions reduced the MD incidence and increased the survival of challenged birds. Furthermore, using the best-fit pathotyping method to rank the virulence, the modified vv+MDV 686 viruses resulted in a pathotype similar to the vvMDV Md5 strain. Thus, these results support our hypothesis that SNPs in one or more of these ORFs are associated with virulence but, as a group, are not sufficient to result in a vMDV pathotype, suggesting that there are additional variants in the MDV genome associated with virulence, which is not surprising given this complex phenotype and our previous finding of additional variants and SNPs associated with virulence.

A pangenome graph reference of 30 chicken genomes allows genotyping of large and complex structural variants.

BACKGROUND: The red junglefowl, the wild outgroup of domestic chickens, has historically served as a reference for genomic studies of domestic chickens. These studies have provided insight into the etiology of traits of commercial importance. However, the use of a single reference genome does not capture diversity present among modern breeds, many of which have accumulated molecular changes due to drift and selection. While reference-based resequencing is well-suited to cataloging simple variants such as single-nucleotide changes and short insertions and deletions, it is mostly inadequate to discover more complex structural variation in the genome. METHODS: We present a pangenome for the domestic chicken consisting of thirty assemblies of chickens from different breeds and research lines. RESULTS: We demonstrate how this pangenome can be used to catalog structural variants present in modern breeds and untangle complex nested variation. We show that alignment of short reads from 100 diverse wild and domestic chickens to this pangenome reduces reference bias by 38%, which affects downstream genotyping results. This approach also allows for the accurate genotyping of a large and complex pair of structural variants at the K feathering locus using short reads, which would not be possible using a linear reference. CONCLUSIONS: We expect that this new paradigm of genomic reference will allow better pinpointing of exact mutations responsible for specific phenotypes, which will in turn be necessary for breeding chickens that meet new sustainability criteria and are resilient to quickly evolving pathogen threats.

An atlas of regulatory elements in chicken: A resource for chicken genetics and genomics.

A comprehensive characterization of regulatory elements in the chicken genome across tissues will have substantial impacts on both fundamental and applied research. Here, we systematically identified and characterized regulatory elements in the chicken genome by integrating 377 genome-wide sequencing datasets from 23 adult tissues. In total, we annotated 1.57 million regulatory elements, representing 15 distinct chromatin states, and predicted about 1.2 million enhancer-gene pairs and 7662 super-enhancers. This functional annotation of the chicken genome should have wide utility on identifying regulatory elements accounting for gene regulation underlying domestication, selection, and complex trait regulation, which we explored. In short, this comprehensive atlas of regulatory elements provides the scientific community with a valuable resource for chicken genetics and genomics.

The immune cell landscape and response of Marek's disease resistant and susceptible chickens infected with Marek's disease virus.

Genetically resistant or susceptible chickens to Marek's disease (MD) have been widely used models to identify the molecular determinants of these phenotypes. However, these prior studies lacked the basic identification and understanding of immune cell types that could be translated toward improved MD control. To gain insights into specific immune cell types and their responses to Marek's disease virus (MDV) infection, we used single-cell RNA sequencing (scRNAseq) on splenic cells from MD resistant and susceptible birds. In total, 14,378 cells formed clusters that identified various immune cell types. Lymphocytes, specifically T cell subtypes, were the most abundant with significant proportional changes in some subtypes upon infection. The largest number of differentially expressed genes (DEG) response was seen in granulocytes, while macrophage DEGs differed in directionality by subtype and line. Among the most DEG in almost all immune cell types were granzyme and granulysin, both associated with cell-perforating processes. Protein interactive network analyses revealed multiple overlapping canonical pathways within both lymphoid and myeloid cell lineages. This initial estimation of the chicken immune cell type landscape and its accompanying response will greatly aid efforts in identifying specific cell types and improving our knowledge of host response to viral infection.

Contribution of the TCRß Repertoire to Marek's Disease Genetic Resistance in the Chicken.

Marek's disease (MD) is a lymphoproliferative disease of chickens induced by Marek's disease virus (MDV), an oncogenic α-herpesvirus. MDV has increased in virulence, prompting continued efforts in both improved vaccines and enhanced genetic resistance. Model pairs of genetically MD-resistant and MD-susceptible chickens that were either MHC-matched or MHC-congenic allowed characterization of T cell receptor (TCR) repertoires associated with MDV infection. MD-resistant chickens showed higher usage of Vß-1 TCRs than susceptible chickens in both the CD8 and CD4 subsets in the MHC-matched model, and in the CD8 subset only in the MHC-congenic model, with a shift towards Vß-1+ CD8 cells during MDV infection. Long and short read sequencing identified divergent TCRß loci between MHC-matched MD-resistant and MD-susceptible chickens, with MD-resistant chickens having more TCR Vß1 genes. TCR Vß1 CDR1 haplotype usage in MD-resistant x MD-susceptible F1 birds by RNAseq indicated that the most commonly used CDR1 variant was unique to the MD-susceptible line, suggesting that selection for MD resistance in the MHC-matched model optimized the TCR repertoire away from dominant recognition of one or more B2 haplotype MHC molecules. Finally, TCR downregulation during MDV infection in the MHC-matched model was strongest in the MD-susceptible line, and MDV reactivation downregulated TCR expression in a tumor cell line.

Fast-forwarding evolution-Accelerated adaptation in a proofreading-deficient hypermutator herpesvirus.

Evolution relies on the availability of genetic diversity for fitness-based selection. However, most deoxyribonucleic acid (DNA) viruses employ DNA polymerases (Pol) capable of exonucleolytic proofreading to limit mutation rates during DNA replication. The relative genetic stability produced by high-fidelity genome replication can make studying DNA virus adaptation and evolution an intensive endeavor, especially in slowly replicating viruses. Here, we present a proofreading-impaired Pol mutant (Y547S) of Marek's disease virus that exhibits a hypermutator phenotype while maintaining unimpaired growth in vitro and wild-type (WT)-like pathogenicity in vivo. At the same time, mutation frequencies observed in Y547S virus populations are 2-5-fold higher compared to the parental WT virus. We find that Y547S adapts faster to growth in originally non-permissive cells, evades pressure conferred by antiviral inhibitors more efficiently, and is more easily attenuated by serial passage in cultured cells compared to WT. Our results suggest that hypermutator viruses can serve as a tool to accelerate evolutionary processes and help identify key genetic changes required for adaptation to novel host cells and resistance to antiviral therapy. Similarly, the rapid attenuation achieved through adaptation of hypermutators to growth in cell culture enables identification of genetic changes underlying attenuation and virulence, knowledge that could practically exploited, e.g. in the rational design of vaccines.

Prediction of transcript isoforms in 19 chicken tissues by Oxford Nanopore long-read sequencing.

To identify and annotate transcript isoforms in the chicken genome, we generated Nanopore long-read sequencing data from 68 samples that encompassed 19 diverse tissues collected from experimental adult male and female White Leghorn chickens. More than 23.8 million reads with mean read length of 790 bases and average quality of 18.2 were generated. The annotation and subsequent filtering resulted in the identification of 55,382 transcripts at 40,547 loci with mean length of 1,700 bases. We predicted 30,967 coding transcripts at 19,461 loci, and 16,495 lncRNA transcripts at 15,512 loci. Compared to existing reference annotations, we found ∼52% of annotated transcripts could be partially or fully matched while ∼47% were novel. Seventy percent of novel transcripts were potentially transcribed from lncRNA loci. Based on our annotation, we quantified transcript expression across tissues and found two brain tissues (i.e., cerebellum and cortex) expressed the highest number of transcripts and loci. Furthermore, ∼22% of the transcripts displayed tissue specificity with the reproductive tissues (i.e., testis and ovary) exhibiting the most tissue-specific transcripts. Despite our wide sampling, ∼20% of Ensembl reference loci were not detected. This suggests that deeper sequencing and additional samples that include different breeds, cell types, developmental stages, and physiological conditions, are needed to fully annotate the chicken genome. The application of Nanopore sequencing in this study demonstrates the usefulness of long-read data in discovering additional novel loci (e.g., lncRNA loci) and resolving complex transcripts (e.g., the longest transcript for the TTN locus).

Identification and Validation of Ikaros (IKZF1) as a Cancer Driver Gene for Marek's Disease Virus-Induced Lymphomas.

Marek's disease virus (MDV) is the causative agent for Marek's disease (MD), which is characterized by T-cell lymphomas in chickens. While the viral Meq oncogene is necessary for transformation, it is insufficient, as not every bird infected with virulent MDV goes on to develop a gross tumor. Thus, we postulated that the chicken genome contains cancer driver genes; i.e., ones with somatic mutations that promote tumors, as is the case for most human cancers. To test this hypothesis, MD tumors and matching control tissues were sequenced. Using a custom bioinformatics pipeline, 9 of the 22 tumors analyzed contained one or more somatic mutation in Ikaros (IKFZ1), a transcription factor that acts as the master regulator of lymphocyte development. The mutations found were in key Zn-finger DNA-binding domains that also commonly occur in human cancers such as B-cell acute lymphoblastic leukemia (B-ALL). To validate that IKFZ1 was a cancer driver gene, recombinant MDVs that expressed either wild-type or a mutated Ikaros allele were used to infect chickens. As predicted, birds infected with MDV expressing the mutant Ikaros allele had high tumor incidences (~90%), while there were only a few minute tumors (~12%) produced in birds infected with the virus expressing wild-type Ikaros. Thus, in addition to Meq, key somatic mutations in Ikaros or other potential cancer driver genes in the chicken genome are necessary for MDV to induce lymphomas.

Marek's Disease Virus Telomeric Integration Profiles of Neoplastic Host Tissues Reveal Unbiased Chromosomal Selection and Loss of Cellular Diversity during Tumorigenesis.

The avian α-herpesvirus known as Marek's disease virus (MDV) linearly integrates its genomic DNA into host telomeres during infection. The resulting disease, Marek's disease (MD), is characterized by virally-induced lymphomas with high mortality. The temporal dynamics of MDV-positive (MDV+) transformed cells and expansion of MD lymphomas remain targets for further understanding. It also remains to be determined whether specific host chromosomal sites of MDV telomere integration confer an advantage to MDV-transformed cells during tumorigenesis. We applied MDV-specific fluorescence in situ hybridization (MDV FISH) to investigate virus-host cytogenomic interactions within and among a total of 37 gonad lymphomas and neoplastic splenic samples in birds infected with virulent MDV. We also determined single-cell, chromosome-specific MDV integration profiles within and among transformed tissue samples, including multiple samples from the same bird. Most mitotically-dividing cells within neoplastic samples had the cytogenomic phenotype of 'MDV telomere-integrated only', and tissue-specific, temporal changes in phenotype frequencies were detected. Transformed cell populations composing gonad lymphomas exhibited significantly lower diversity, in terms of heterogeneity of MDV integration profiles, at the latest stages of tumorigenesis (>50 days post-infection (dpi)). We further report high interindividual and lower intraindividual variation in MDV integration profiles of lymphoma cells. There was no evidence of integration hotspots into a specific host chromosome(s). Collectively, our data suggests that very few transformed MDV+ T cell populations present earlier in MDV-induced lymphomas (32-50 dpi), survive, and expand to become the dominant clonal population in more advanced MD lymphomas (51-62 dpi) and establish metastatic lymphomas.

3D genomics across the tree of life reveals condensin II as a determinant of architecture type.

We investigated genome folding across the eukaryotic tree of life. We find two types of three-dimensional (3D) genome architectures at the chromosome scale. Each type appears and disappears repeatedly during eukaryotic evolution. The type of genome architecture that an organism exhibits correlates with the absence of condensin II subunits. Moreover, condensin II depletion converts the architecture of the human genome to a state resembling that seen in organisms such as fungi or mosquitoes. In this state, centromeres cluster together at nucleoli, and heterochromatin domains merge. We propose a physical model in which lengthwise compaction of chromosomes by condensin II during mitosis determines chromosome-scale genome architecture, with effects that are retained during the subsequent interphase. This mechanism likely has been conserved since the last common ancestor of all eukaryotes.

Functional annotations of three domestic animal genomes provide vital resources for comparative and agricultural research.

Gene regulatory elements are central drivers of phenotypic variation and thus of critical importance towards understanding the genetics of complex traits. The Functional Annotation of Animal Genomes consortium was formed to collaboratively annotate the functional elements in animal genomes, starting with domesticated animals. Here we present an expansive collection of datasets from eight diverse tissues in three important agricultural species: chicken (Gallus gallus), pig (Sus scrofa), and cattle (Bos taurus). Comparative analysis of these datasets and those from the human and mouse Encyclopedia of DNA Elements projects reveal that a core set of regulatory elements are functionally conserved independent of divergence between species, and that tissue-specific transcription factor occupancy at regulatory elements and their predicted target genes are also conserved. These datasets represent a unique opportunity for the emerging field of comparative epigenomics, as well as the agricultural research community, including species that are globally important food resources.

Depletion of CD8αß+ T Cells in Chickens Demonstrates Their Involvement in Protective Immunity towards Marek's Disease with Respect to Tumor Incidence and Vaccinal Protection.

Marek's disease (MD) is a lymphoproliferative disease in chickens caused by Marek's disease virus (MDV), a highly oncogenic alphaherpesvirus. Since 1970, MD has been controlled through widespread vaccination of commercial flocks. However, repeated and unpredictable MD outbreaks continue to occur in vaccinated flocks, indicating the need for a better understanding of MDV pathogenesis to guide improved or alternative control measures. As MDV is an intracellular pathogen that infects and transforms CD4+ T cells, the host cell-mediated immune response is considered to be vital for controlling MDV replication and tumor formation. In this study, we addressed the role of CD8+ T cells in vaccinal protection by widely-used monovalent (SB-1 and HVT) and bivalent (SB-1+HVT) MD vaccines. We established a method to deplete CD8+ T cells in chickens and found that their depletion through injection of anti-CD8 monoclonal antibodies (mAb) increased tumor induction and MD pathology, and reduced vaccinal protection to MD, which supports the important role of CD8+ T cells for both MD and vaccinal protection.

Genome-wide characterization of copy number variations in the host genome in genetic resistance to Marek's disease using next generation sequencing.

BACKGROUND: Marek's disease (MD) is a highly neoplastic disease primarily affecting chickens, and remains as a chronic infectious disease that threatens the poultry industry. Copy number variation (CNV) has been examined in many species and is recognized as a major source of genetic variation that directly contributes to phenotypic variation such as resistance to infectious diseases. Two highly inbred chicken lines, 63 (MD-resistant) and 72 (MD-susceptible), as well as their F1 generation and six recombinant congenic strains (RCSs) with varied susceptibility to MD, are considered as ideal models to identify the complex mechanisms of genetic and molecular resistance to MD. RESULTS: In the present study, to unravel the potential genetic mechanisms underlying resistance to MD, we performed a genome-wide CNV detection using next generation sequencing on the inbred chicken lines with the assistance of CNVnator. As a result, a total of 1649 CNV regions (CNVRs) were successfully identified after merging all the nine datasets, of which 90 CNVRs were overlapped across all the chicken lines. Within these shared regions, 1360 harbored genes were identified. In addition, 55 and 44 CNVRs with 62 and 57 harbored genes were specifically identified in line 63 and 72, respectively. Bioinformatics analysis showed that the nearby genes were significantly enriched in 36 GO terms and 6 KEGG pathways including JAK/STAT signaling pathway. Ten CNVRs (nine deletions and one duplication) involved in 10 disease-related genes were selected for validation by using quantitative real-time PCR (qPCR), all of which were successfully confirmed. Finally, qPCR was also used to validate two deletion events in line 72 that were definitely normal in line 63. One high-confidence gene, IRF2 was identified as the most promising candidate gene underlying resistance and susceptibility to MD in view of its function and overlaps with data from previous study. CONCLUSIONS: Our findings provide valuable insights for understanding the genetic mechanism of resistance to MD and the identified gene and pathway could be considered as the subject of further functional characterization.

Pathogen transmission from vaccinated hosts can cause dose-dependent reduction in virulence.

Many livestock and human vaccines are leaky because they block symptoms but do not prevent infection or onward transmission. This leakiness is concerning because it increases vaccination coverage required to prevent disease spread and can promote evolution of increased pathogen virulence. Despite leakiness, vaccination may reduce pathogen load, affecting disease transmission dynamics. However, the impacts on post-transmission disease development and infectiousness in contact individuals are unknown. Here, we use transmission experiments involving Marek disease virus (MDV) in chickens to show that vaccination with a leaky vaccine substantially reduces viral load in both vaccinated individuals and unvaccinated contact individuals they infect. Consequently, contact birds are less likely to develop disease symptoms or die, show less severe symptoms, and shed less infectious virus themselves, when infected by vaccinated birds. These results highlight that even partial vaccination with a leaky vaccine can have unforeseen positive consequences in controlling the spread and symptoms of disease.

Towards a mechanistic understanding of the synergistic response induced by bivalent Marek's disease vaccines to prevent lymphomas.

BACKGROUND: Marek's disease (MD) is a lymphoproliferative disease of chickens caused by Marek's disease virus (MDV), an oncogenic α-herpesvirus. Since 1970, MD has been controlled by widespread vaccination; however, more effective MD vaccines are needed to counter more virulent MDV strains. The bivalent vaccine combination of SB-1 and herpesvirus of turkey (HVT) strain FC126 has been widely used. Nonetheless, the mechanism(s) underlying this synergistic effect has not been investigated. METHODS: Three experiments were conducted where SB-1 or HVT were administered as monovalent or bivalent vaccines to newly hatched chickens, then challenged five days later with MDV. In Experiment 1, levels of MDV replication in PBMCs were measured over time, and tumor incidence and vaccinal protection determined. In Experiment 2, MDV and vaccine strains replication levels in lymphoid organs were measured at 1, 5, 10, and 14 days post-challenge (DPC). In Experiment 3, to verify that the bursa was necessary for HVT protection, a subset of chicks were bursectomized and these birds plus controls were similarly vaccinated and challenged, and the levels of vaccinal protection determined. RESULTS: The efficacy of bivalent SB-1 + HVT surpasses that of either SB-1 or HVT monovalent vaccines in controlling the level of pathogenic MDV in PBMCs until the end of the study, and this correlated with the ability to inhibit tumor formation. SB-1 replication in the spleen increased from 1 to 14 DPC, while HVT replicated only in the bursa at 1 DPC. The bursa was necessary for immune protection induced by HVT vaccine. CONCLUSION: Synergy of SB-1 and HVT vaccines is due to additive influences of the individual vaccines acting at different times and target organs. And the bursa is vital for HVT to replicate and induce immune protection.

Identification of Marek's disease virus genes associated with virulence of US strains.

Marek's disease virus (MDV) is the most well-cited example of vaccine-driven virulence evolution. MDV induces a lymphoproliferative disease in chickens, which is currently controlled by widespread vaccination of flocks. Unfortunately, Marek's disease (MD) vaccines, while effective in preventing tumours, do not prevent viral replication and mutation, which has been hypothesized as the major driving force for increased MDV virulence of field strains during the past 40 years in US commercial flocks. To limit future virulence increases, there is interest in characterizing MDV strain genomes collected over the years and associating genetic variations with variation in virulence. In this study, we characterized 70 MDV genomes with known virulence by complete or targeted DNA sequencing, and identified genetic variants that showed association with virulence. Our results revealed a number of MDV genes as would be expected for a complex trait. In addition, phylogenetic analysis revealed a clear separation of strains that varied by virulence. Interestingly, high virulence isolates from the same farms persisted over years despite eradication attempts, which has implications on control efforts. Given the growing ability to bioengineer the MDV genome, it should be feasible to experimentally test whether these individual variants influence virulence markers alone or combinations. Once validated, these markers may provide an alternative to live bird testing for evaluating virulence of new MDV field strains.

Vaccination and Host Marek's Disease-Resistance Genotype Significantly Reduce Oncogenic Gallid alphaherpesvirus 2 Telomere Integration in Host Birds.

Marek's disease (MD) is an infectious disease characterized by lymphomas and high mortality in susceptible chickens. The causative and ubiquitous alpha-herpesvirus known as MD virus (MDV) integrates into host telomeres during early infection through latency, known to be an important phase for oncogenic transformation. Herein, we sought to determine the influence of vaccination and host genetics on the temporal dynamics of MDV-host genome interactions. We studied integration profiles using 2 MD vaccines that vary in protective efficacy in 2 genetic lines that differ in MD resistance/susceptibility. Virus integration of both oncogenic MDV and vaccine strains was observed in both MD susceptible and resistant birds, however, the lines differed in their dynamic telomere-integration profiles. Notably, the resistant host genotype exhibited a smaller percentage of replicating cells with the virus telomere-integrated only phenotype as compared to the susceptible genotype. Vaccination with Rispens, the most protective MD vaccine, also reduced the establishment of the virus telomere-integrated only phenotype, suggesting a significant role of the phenotype in MD lymphoma development. The effect of Rispens vaccination was most dramatic in the susceptible genotype. These results suggest important connections between vaccinal immunity, MDV telomere integration, virus-induced oncogenesis, and virus-host genome interactions in the context of host genetics and disease susceptibility.

Genetic assessment of inbred chicken lines indicates genomic signatures of resistance to Marek's disease.

BACKGROUND: Marek's disease (MD) is a highly contagious pathogenic and oncogenic disease primarily affecting chickens. However, the mechanisms of genetic resistance for MD are complex and not fully understood. MD-resistant line 63 and MD-susceptible line 72 are two highly inbred progenitor lines of White Leghorn. Recombinant Congenic Strains (RCS) were developed from these two lines, which show varied susceptibility to MD. RESULTS: We investigated genetic structure and genomic signatures across the genome, including the line 63 and line 72, six RCSs, and two reciprocally crossed flocks between the lines 63 and 72 (F1 63 × 72 and F1 72 × 63) using Affymetrix® Axiom® HD 600 K genotyping array. We observed 18 chickens from RCS lines were specifically clustered into resistance sub-groups distributed around line 63. Additionally, homozygosity analysis was employed to explore potential genetic components related to MD resistance, while runs of homozygosity (ROH) are regions of the genome where the identical haplotypes are inherited from each parent. We found several genes including SIK, SOX1, LIG4, SIK1 and TNFSF13B were contained in ROH region identified in resistant group (line 63 and RCS), and these genes have been reported that are contribute to immunology and survival. Based on FST based population differential analysis, we also identified important genes related to cell death and anti-apoptosis, including AKT1, API5, CDH13, CFDP and USP15, which could be involved in divergent selection during inbreeding process. CONCLUSIONS: Our findings offer valuable insights for understanding the genetic mechanism of resistance to MD and the identified genes could be considered as candidate biomarkers in further evaluation.

Evaluation and Identification of Marek's Disease Virus BAC Clones as Standardized Reagents for Research.

Marek's disease virus (MDV) is an alphaherpesvirus that causes Marek's disease (MD), a lymphoproliferative disease in chickens. Understanding of MDV gene function advanced significantly following the cloning of the MDV genome as either a series of overlapping cosmids or as a bacterial artificial chromosome (BAC), both of which could produce viable MDV. The objectives of this study were to compare multiple virulent MDV BAC clones using the Avian Disease and Oncology Laboratory's pathotyping assay, and to demonstrate the use of these clones as standardized reagents for a modified pathotyping assay by other laboratories. To date, MDV BAC clones have been produced for at least 10 MDV strains from all three serotypes including several virulent serotype 1 strains. We determined that MDV BAC clones exist for each virulent pathotype, despite the fact that these clones are not always equal in virulence to their corresponding parental strains. One clone from each pathotype was further evaluated in commercial specific-pathogen-free (SPF) chickens and found suitable for use in assays such as best-fit pathotyping, although results were variable based on the source of the SPF birds. The benefits of using BAC clones, which include easy shipping, ability to more easily manipulate, and long-term ability to use at a low passage level, are likely to result in the use of BAC clones as standard reagents for MD research. The use of the defined set of clones should allow side-by-side comparison, allowing researchers to better interpret results produced in different laboratories using different MDV field strains. Furthermore, a modified best-fit pathotyping assay has been proposed using these clones and reduced bird numbers.